Spin-Orbital Liquid on a Triangular Lattice 
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Using Lanczos exact diagonalization of finite clusters we demonstrate that the spin-orbital d 
model for triply degenerate ti g orbitals on a triangular lattice provides an example of a spin-orbital 
liquid ground state. We also show that the spin-orbital liquid involves entangled valence bond states 
which violate the Goodenough-Kanamori rules, and modify effective spin exchange constants. 
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I. INTRODUCTION 

A variety of very interesting and challenging problems 
in condensed matter theory arises in systems of strongly 
correlated electrons with degenerate orbitals When 
the intraorbital Coulomb element U is much larger than 
the effective electron hopping i, i.e., U t, the magnetic 
properties follow from the spin-orbital superexchange Q ■ 
Usually spin interactions are then determined by orbital 
order and complementary types of spin and orbital or- 
der coexist in agreement with the Goodenough-Kanamori 
rules (GKR) [3]. However, large quantum fluctuations 
that emerge from strongly frustrated orbital interactions 
could instead stabilize disordered phases [4j. This ob- 
servation triggered the search for an example of a spin- 
orbital liquid (SOL) ground state (GS), similar to a spin 
liquid state in a one-dimensional Heisenberg antiferro- 
magnet. In principle such a SOL state might be expected 
for the spin-orbital d 9 model on a triangular lattice in 
LiNi02 , but Ising-like orbital interactions suppress it [j| . 

Exotic behavior of spin-orbital systems follows from 
spin-orbital entanglement (SOE) Q . To name a few phe- 
nomena, it is responsible for the temperature dependence 
of optical spectral weights in LaV03 [7j , plays a role in 
the phase diagram of the -RVO3 perovskites H , and re- 
stricts propagation of a hole in states with entangled spin- 
orbital order Recently novel phases with SOE were 
discovered in a bilayer spin-orbital d 9 model [lol ] , but also 
here a SOL phase could not be established. 



II. MODEL AND CALCULATION METHOD 



The orbital interactions for t2 g orbitals, with T =1/2 
pseudospins that depend on the bond direction, are more 
quantum than for e g ones and all three pseudospin com- 
ponents contribute for each bond [if. Here we em- 
ploy Lanczos exact diagonalization to investigate a spin- 
orbital d 1 model for triply degenerate ti g orbitals on a 
triangular lattice as in NaTi02 [H|, with superexchange 
(jhLs)i direct exchange (Hd) and mixed terms responsible 



for enhanced quantum fluctuations (H m ), 

U = j[{l-a)U s + y/(l - a)a U m + a H d ] . (1) 

The model depends on the parameter < a < 1 and 
on Hund's exchange 77 = Jjj/U. Here J is the exchange 
energy. In the direct exchange limit (a — I) the model 
Eq. (1) is exactly solvable and the GS was determined 
by considering the valence bond (VB) dimcr coverings 
of the lattice with each dimer containing a spin singlet 
[l2| . Below we show by analyzing the results of Lanczos 
diagonalization 

.13] that a SOL is realized in the opposite 
superexchange limit (a — 0) of the model. 

The essential feature of the model Eq. (fTJ) is that 
5=1/2 quantum spins are coupled by an SU(2) sym- 
metric (Si ■ Sj) interactions, while the orbital interac- 
tions obey a much lower discrete symmetry between three 
equivalent directions 7 = a, b, c in the lattice, with two tig 

orbital flavors active in the superexchange via (Tj -Tj)^ 
term for T = 1/2 pseudospin, and the third one con- 
tributing to direct exchange; for more details see [ll|, [H[ . 
The drastic difference in occupied orbital states realized 
in both limits at 77 = is illustrated for the case of a 
9-site cluster by insets in Fig. 1(b) — one finds equal 
occupancies of each orbital state at a = 0, while 756 de- 
generate orbital dimer VB-like states are found at a = 1 
[I2I [l3l]. and the degeneracy scales with system size. 

We characterize the GS by spin, orbital and spin- 
orbital (four-operator) bond correlations (d is the degen- 
eracy of the GS {|n)}), given by 



S = \ ZK n l^'^'l n ) 

n 



(2) 
(3) 



±. y £(n\0 i .§ j -S)(? i .f j -T)™\n). (4) 



Below we consider two clusters with periodic boundary 
conditions (PBC): a hexagonal cluster of N = 7 sites 



2 



0.0 








I 














( 








a — & e 


-0.1 




























-0.2 
-0.3 


s 






T 


h-B 


—I 


C 



















(b) 

0.0 
<o -0.1 
<0 -0.2 
-0.3 









? 1 1 



0.0 0.2 0.4 0.6 0.8 1.0 

a 



FIG. 1: Bond correlations at 77 = 0: spin (S), orbital (T), and 
spin-orbital (C) for: (a) N7 cluster, and (b) N9 cluster. The 
vertical lines indicate the exactly determined range of possible 
values that follows from the GS degeneracy. The insets in (b) 
indicate typical orbital patterns in the superexchange (a = 0) 
and direct exchange (a = 1) limit for the rhombic N9 cluster. 



(N7) and a rhombic cluster of N = 9 sites (N9). Due 
to PBC all sites and bonds are equivalent and the above 
correlations do not depend on the bond (ij) and its di- 
rection 7. Each tig orbital is occupied on average by 1/3 
electron, but the states are manifestly different in the 
limits of a = and a = 1, see the insets in Fig. 1(b). 



III. RESULTS AND DISCUSSION 

In both N7 and N9 cluster spin S and orbital T cor- 
relations are negative and the GKR stating that these 
correlations should be complementary are violated, see 
Fig. 1. Frustration in the triangular lattice is responsible 
for a rather weak and independent of a antiferromagnetic 
(AF) spin correlations in the N7 cluster, S ~ —0.107, see 
Fig. 1(a). These correlations are weaker (S ~ —0.090) 
at a — in the N9 cluster and become more pronounced 
(S ~ —0.144) when a ~ 0.6 and joint spin-orbital fluc- 
tuations weaken to C ~ —0.050, see Fig. 1(b). The or- 
bital correlations weaken as well when a increases toward 
a = 1, particularly in the N9 cluster. Joint spin-orbital 
correlations are similar in both clusters (e.g. C ~ —0.070 
at a. = 0) and \C\ gradually decreases when spin and 
orbitals disentangle approaching a = 1. 

An important question is whether spin order and ex- 
citations could be described by an effective spin model 



derived from Eq. ([]}. In order to illustrate consequences 
of SOE in magnetic states we rewrite the d 1 spin-orbital 
model Eq. (|TJ) in a general form [2j resembling a spin 
Hamiltonian, 

n = E {4 7) (^-^)+^ 7) }' ( 5 ) 

where the orbital operators and depend on the 
parameters {a, 77} for a bond (ij) along axis 7. Mean 
field (MF) procedure used frequently reads [l3j], 

+ E {4 7) • ^> + ■ 

<*j) IIt 

It disentangles spin and orbital degrees of freedom and 
is used to determine the MF spin constant for N7 and 
N9 clusters by averaging the orbital operator (its 
explicit form is given in [13j) over the MF GS l^o), 

J MF = ($ol4 7) l*o)- ( ? ) 

Note that the orbital fluctuations in the term oc 
\/a(l — a) in Eq. (1) contribute here as well as they 
couple different components of |$o)- I n contrast, the ex- 
act exchange constant, J eX act, is found when the exact 
GS |$) obtained after Lanczos diagonalization is used. 

In Fig. 2 we compare the phase diagrams obtained 
from the above MF procedure and from exact diagonal- 
ization for the N9 cluster. Consider first a quantum phase 
transition from the low-spin (St — 1/2) disordered phase 
to the high-spin (St — 9/2) ferromagnetic (FM) phase 
which occurs for sufficiently large 77. When spin and or- 
bital operators are disentangled in Eq. ©, i.e., C = 
Q , it coincides with the sign change of the MF exchange 
constant Jmf and no other phase (with l/2<5t<9/2) 
is found, see Fig. 2(a), as in a spin system. 

Comparing the values of Jmf and J eX act found from 
the MF and from exact diagonalization of the N9 cluster 
(Fig. 2), one finds that J exa ct > Jmf in a broad regime 
of a except for a ~ 1. Therefore, the MF approximation 
turns out to be rather unrealistic and overestimates (un- 
derestimates) the stability of states with FM (AF) spin 
correlations. The value of Jmf decreases with increas- 
ing 77, but positive values Jmf > are found at 77 = 
only if 0.07 < a < 1. This demonstrates that FM states: 
(i) are favored when joint spin-orbital fluctuations are 
suppressed, and (ii) are stabilized by orbital fluctuations 
close to a = even in absence of Hund's exchange. The 
transition from the low-spin (St = 1/2) to the high-spin 
(Sf = 9/2) state occurs in presence of SOE at a much 
higher value of 77 rs 0.14, with only weak dependence on 
a, see Fig. 2b. In addition, one finds a phase with an 
intermediate spin value S t — 3/2 for < a < 0.21 and 
0.44 < a < 0.88, and J eX act changes discontinuously at 
the transition to the FM phase. 
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(a) (1-3t|)Jmf/J 




(b) (l-3r))Je Xact /J 
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FIG. 2: Phase diagram in the (a, n) plane (heavy lines) and 
exchange constants (contour plots) as obtained for the N9 
cluster with PBC: (a) within the MF calculation, and (b) 
using exact Lanczos diagonalization. In the MF case (a) the 
transition from low-spin (St = 1/2) to high-spin (<S t = 9/2) 
phase occurs when Jmf changes sign. In the exact calculation 
one finds in addition an intermediate phase with St = 3/2, 
stable between the heavy black lines in (b). 

Altogether, the qualitative trends found for the N9 
cluster are generic and agree with those observed for the 
N7 cluster, see Fig. 15 in [L3J. In both cases one finds 
that: (i) the FM phase is stable in the MF approxima- 
tion close to a = and becomes degenerate with the 
low-spin phase at a = 1, (ii) the MF procedure is exact 
in the regime of FM phase, and (Hi) the transition to the 
FM phase occurs gradually through intermediate values 
of total spin St (except at a = 1). This suggests that 
partially polarized FM phase should occur in the ther- 
modynamic limit of the model Eq. (1) and arises due to 
SOE which is gradually suppressed when rj increases. 

We argue that the present study and the results pre- 
sented in [l3[ provide evidence in favor of a quantum 



SOL phase in the present d 1 spin-orbital model Eq. (1) 
in a corner of its phase diagram — in the regime of small 
values of both a and 77 parameter. In agreement with the 
directional nature of orbital interactions, this SOL phase 
is also characterized by rather strong VB dimer corre- 
lations 14J. The consequences of SOE are more severe 
in this phase and the transition to the FM phase occurs 
typically at a much higher value of r\ than the one where 
•/exact changes its sign. Therefore, we suggest that even 
in case when magnetic exchange J oxa ct is accurately eval- 
uated using the relevant orbital correlations, it looses its 
predictive power and is inadequate to describe the mag- 
netic ground state and excitations in a system dominated 
by SOE, where the GKR are violated. Note that the frus- 
trated triangular lattice plays here an important role and 
removes any kind of orbital order. 



IV. SUMMARY AND CONCLUSIONS 

Summarizing, we have demonstrated that the GKR are 
violated in the regime of weak Hund's exchange and the 
spin-orbital liquid phase is stabilized by spin-orbital en- 
tanglement in the d 1 spin-orbital model on the triangular 
lattice. In this regime the MF decoupling procedure of 
spin and orbital operators fails and the magnetic proper- 
ties can be determined only by solving the full entangled 
spin-orbital many-body problem. 

Finally, we point out that spin-orbital entangled states 
play a role in exotic ground states [l(| as well as low en- 
ergy excitations for realistic values of r\ ~ 0.14. Their 
consequences have already been established in the vana- 
dium perovskites 043 , and we expect that they could be 
of even more importance in systems with geometrically 
frustrated lattice such as the one considered here. 
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